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Abstract. Dimuon production has been studied in a series of fixed-target experiments at Fermilab 
during the last two decades. Highlights from these experiments, together with recent results from 
■ the Fermilab E866 experiment, are presented. Future prospects for studying the parton distributions 

^ ' in the nucleons and nuclei using dimuon production are also discussed. 
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INTRODUCTION 



The Drell-Yan process [1], in which a charged lepton pair is produced in a hadron-hadron 
interaction via the electromagnetic qq — )■ /+/^ process, has provided unique information 
^1 ■ on parton distributions. In particular, the Drell-Yan process has been used to determine 

the antiquark contents of nucleons and nuclei [2], as well as the quark distributions of 
pions, kaons, and antiprotons [3]. Such information is difficult, if not impossible, to 
^ ■ obtain from DIS experiments. As the Drell-Yan process can be well described by next- 

^ . to-leading order QCD calculations [4], a firm theoretical framework exists for utilizing 

O ! the Drell-Yan process to extract the parton distributions. 

■ During the last two decades, a series of fixed-target dimuon production experiments 

. (E772, E789, E866) have been carried out using 800 GeV/c proton beam at Fermilab. 

At 800 GeV/c, the dimuon data contain Drell-Yan continuum up to dimuon mass of 
~ 15 GeV as well as quarkonium productions (J/^, and T resonances). The Drell- 
Yan process and quarkonium productions often provide complementary information, 
since Drell-Yan is an electromagnetic process via quark-antiquark annihilation while 
the quarkonium production is a strong interaction process dominated by gluon-gluon 
fusion at this beam energy. 

The Fermilab dimuon experiments covers a broad range of physics topics. The Drell- 
Yan data have provided informations on the antiquark distributions in the nucleons [5, 
6, 7, 8] and nuclei [9, 10]. These results showed the surprising results that the antiquark 
distributions in the nuclei are not enhanced [9, 10], contrary to the predictions of models 
which explain the EMC effect in term of nuclear enhancement of exchanged mesons. 
Moreover, the Drell-Yan cross section ratios p + d/ p + p clearly establish the flavor 
asymmetry of the d and u distributions in the proton, and they map out the x-dependence 
of this asymmetry [6, 7, 8]. Pronounced nuclear dependences of quarkonium productions 
have been observed for J/^, and T resonances [11, 12, 13, 14]. 

Several review articles covering some of these results are available [2, 15, 16]. In this 
article, we will focus on the recent results from experiment E866 and future prospect of 
dimuon experiments at Fermilab and J-PARC. 




FLAVOR STRUCTURE OF LIGHT-QUARK SEA 

In the CERN NA51 [17] and Fermilab E866 [6, 7, 8] experiments on proton-induced 
dimuon production, a striking difference was observed for the Drell-Yan cross sec- 
tions between p + p and p + d. As the underlying mechanism for the Drell-Yan pro- 
cess involves quark-antiquark annihilation, this difference has been attributed to the 
asymmetry between the up and down sea quark distributions in the proton. From the 
o{p + d)i)Y /'2.o{p + p)dy ratios the Bjorken-x dependence of the sea-quark d/u flavor 
asymmetry has been extracted [6, 7, 8, 17]. 

Future fixed-target dimuon experiments have been proposed at the 120 GeV Fermilab 
Main Injector (FMI) and the 50 GeV J-PARC facilities. The Fermilab proposal [18], 
E906, has been approved and is expected to start data-taking around 201 1. Two dimuon 
proposals (P04 [19] and P24 [20]) have also been submitted to the J-PARC for approval. 
The lower beam energies at FMI and J-PARC present opportunities for extending the 
d/u and the nuclear antiquark distribution measurements to larger x {x > 0.25). For 
given values of x\ and X2, the Drell-Yan cross section is proportional to l/s, hence 
a gain of ~ 16 times in the Drell-Yan cross sections can be obtained at the J-PARC 
energy of 50 GeV. Since the perturbative process gives a symmetric d/u while non- 
perturbative processes are necessary to generate an asymmetric d/u sea, it would be 
very important to extend the Drell-Yan measurements to kinematic regimes beyond the 
current limits. Another advantage of lower beam energies is that a much more sensitive 
study of the partonic energy loss in nuclei could be carried out using the Drell-Yan 
nuclear dependence [21]. 
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FIGURE 2. The E866 CT{p + d)/2o{p + p) cross section ratios for T resonances as a function of X2. 
The corresponding ratios for Drell- Yan cross sections are also shown. The error bars are statistical only. 



To disentangle the d/u asymmetry from the possible charge-symmetry violation ef- 
fect [22, 23], one could consider W boson production in p-\-p collision at RHIC. An in- 
teresting quantity to be measured is the ratio of the p + p ^ +x and p + p ^ W +x 
cross sections [24]. It can be shown that this ratio is very sensitive to d/u. An important 
feature of the W production asymmetry in p + p collision is that it is completely free 
from the assumption of charge symmetry. Another advantage of production in p + p 
collision is that no nuclear effects need to be considered. Finally, the W production is 
sensitive to d/u flavor asymmetry at a scale of ~ 6500 GeV^/c^, significantly larger 
than all existing measurements. This offers the opportunity to examine the QCD evolu- 
tion of the sea-quark flavor asymmetry. Figure 1 shows the predictions [25] forp+p 
collision at ^/s = 500 GeV. The MRS SO' corresponds to the d/ u symmetric parton dis- 
tributions, while the other three parton distribution functions are from recent global fits 
with asymmetric d/u sea-quark distributions. Figure 1 clearly shows that W asymmetry 
measurements at RHIC could provide an independent determination of d/ u. First results 
from the LHC are also expected soon, and will provide additional information [25]. 

Unlike the electromagnetic Drell- Yan process, quarkonium production is a strong 
interaction dominated by the subprocess of gluon-gluon fusion at 800 GeV beam energy. 
Therefore, the quarkonium production cross sections are primarily sensitive to the gluon 
distributions in the colliding hadrons. The T production ratio, o{p -\- d ^T) /2o{p -\- 
p —^T), is expected to probe the gluon content in the neutron relative to that in the 
proton [26]. 

The o{p + d)/2o{p + p) ratios for T(15 + 25 + 35) production are shown in Fig. 2 
as a function of X2 [27]. These ratios are consistent with unity, in striking contrast to 
the corresponding values for the Drell- Yan process, also shown in Fig. 2. The difference 
between the Drell- Yan and the T cross section ratios clearly reflect the different underly- 
ing mechanisms in these two processes. For T production, the dominance of the gluon- 



gluon fusion subprocess at this beam energy implies that a{p + J — > T) /2a{p + p ^ 
T) pa 5(1 + gni-Xij/gpixj))- Figure 2 shows that the gluon distributions in the proton 
(gp) and neutron (gn) are very similar over the X2 range 0.09 <X2 < 0.25. The overall 
a{p + d ^ T) /2o{p + p ^ T) ratio, integrated over the measured kinematic range, is 
0.984 ± 0.026(stat.) ± 0.01 (syst.) . The Y data indicate that the gluon distributions in the 
proton and neutron are very similar. 

The upcoming Fermilab E906 experiment is expected to provide a precise measure- 
ment of the o{p + d J/^)/2a{p + p J /^) ratio. These data should further test 
the equality of the gluon distributions in the proton and neutron. Moreover, these data 
could identify EMC effects for deuteron. It is interesting to note that analogous ratios 
measured in DIS and the Drell-Yan can not be used to determine the EMC effects for 
deuteron, since different parton distributions are being probed for the hydrogen and the 
deuterium targets. If one assumes that the gluon distributions are identical in proton and 
neutron, then any deviation of the a(p + J — )> //^) /2o{p+ p ^ J /^) ratio from unity 
would be attributed to the EMC effects of gluon distributions in deuteron. 

FLAVOR DEPENDENCE OF THE EMC EFFECT 

Despite a quarter century of significant experimental and theoretical effort, the specific 
origins of the observed A dependence of the nuclear quark distributions have yet to 
be unambiguously identified. Attempts to explain the EMC effect have led to a large 
collection of theoretical models [28, 29], many of which are capable of describing the 
essential features of the data, however the underlying physics mechanisms in each model 
are often very different. 

A new calculation of the modifications of nucleon quark distributions in the nuclear 
medium has recently been reported [30, 31]. In this approach by Cloet, Bentz and 
Thomas (CBT), the Nambu-Jona-Lasinio model is used to describe the coupling of the 
quarks in the bound nucleons to the scalar and vector mean fields inside a nucleus. These 
nucleon quark distributions are then convoluted with a nucleon momentum distribution 
in the nucleus to generate the nuclear quark distributions [30]. An important feature of 
this model is that for N ^ Z nuclei (where A'^ and Z refer to the number of neutrons 
and protons) the i so vector-vector mean field (usually denoted by p^) will affect the up 
quarks differently from the down quarks in the bound nucleons. Therefore, this model 
has a novel prediction that the u and d quarks have distinct nuclear modifications for 
N nuclei. 

Semi-inclusive DIS (SIDIS) on heavy nuclear targets, in which the flavor of the struck 
quark is tagged by the detected hadron, is a promising experimental tool to search for 
the flavor-dependent EMC effect. Recently, Lu and Ma [32] pointed out that charged 
lepton SIDIS off nuclear targets and the deuteron can be used to probe the flavor content 
of the nuclear quark sea, which can help distinguish between the various models of 
the EMC effect. Indeed, a SIDIS experiment [33] aiming at a precise determination of 
flavor dependence of the EMC effect has also been proposed at the upgraded 12GeV 
JLab facility. 

The flavor dependence of the EMC effect is a promising experimental observable 
to distinguish among the plethora of models that can describe the EMC effect. As 





FIGURE 3. The existing data for the ratios 
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discussed earlier, the nuclear dependence of proton-induced Drell-Yan process measured 
at Fermilab has shown the surprising results that the antiquark distributions in the nuclei 
are not enhanced [9, 10], contrary to the predictions of many EMC models. Pion-induced 
Drell-Yan process provides another experimental tool with which search for flavor- 
dependent effects in the nuclear modification of the nucleon structure functions [34]. 
To explore the sensitivity of pion-induced Drell-Yan processes to a flavor-dependent 

EMC effect, we consider the three ratios ^^rj^-^^^j ^ o'^^k-+d) ^"^^ '^'^^^^ ^ 

represents a nuclear, D a deuteron and H a hydrogen target. Assuming isospin symmetry. 



which implies Uj^+ = d^-, Uj^- = dj^+, Uj^+ = d^-, u^- = dj^+ and keeping only the 
dominant terms in each cross-section, one readily obtains 
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The up and down nuclear quark distributions are labeled by ua and respectively, w/j is 
the up quark distribution in the deuteron and Up the up quark distribution in the proton. 
Eqs. 1-3 demonstrate that these Drell-Yan cross-section ratios are very sensitive to the 
flavor dependence of the EMC effect. 

Fig. 3 shows the comparison [34] of the calculations of the pion-induced Drell- 
Yan cross-section ratios with the existing data. The top left panel shows the ratio of 
a^^ in^ +W) 

g^^(;r-+£>) NAIO experiment [35]. These plots contain both the Pbeam = 286 

and 140GeV data sets, which are very similar. The calculations are performed at Pbeam = 
286 GeV, since most of the data was obtained at this energy. The PDFs of the CBT 
model [30, 31] at a fixed of 25 GeV ^, which is approximately the mean of 

the NAIO experiment, were used. The top right panel shows the ratio ^orl^-^^^j from 

the NA3 experiment [36]. The data were collected using a 150 GeV Tt" beam and the 
range covered was 16.8 < < 70.6 GeV^. The calculations are performed for 
Pbeam = 150 GeV and = ISG&V^. 

The solid curves in Fig. 3 are calculations using the flavor-dependent ua{x) and dA{x) 
from the CBT model with N/Z = 1.5, corresponding approximately to the N/Z values 
for the Au, W and Pt nuclei. The dashed curves correspond to the calculated ratios 
using the nuclear PDFs from the CBT model for N = Z. Since ua/uq = dA/do in 
this case, the dashed curves are representative of the predictions for flavor-independent 
EMC models. Figure 3 shows that the NAIO data do not exhibit a clear preference for 
the flavor-dependent versus flavor-independent nuclear PDFs. In contrast, the NA3 data 
strongly favor the calculations using flavor-dependent nuclear PDFs. The existing Drell- 
Yan data are not sufficiently accurate yet, although the NA3 data clearly favor the flavor- 
dependent over the flavor-independent nuclear PDFs. Precise future pion-induced Drell- 
Yan experiments can provide unique constraints that will help distinguish the various 
theoretical models and most importantly shed new light on the origins of the EMC effect. 



TRANSVERSE SPIN AND DRELL-YAN PROCESS 

The study of the transverse momentum dependent (TMD) parton distributions of the 
nucleon has received much attention in recent years as it provides new perspectives on 
the hadron structure and QCD [37]. These novel TMDs can be extracted from semi- 
inclusive deep-inelastic scattering (SIDIS) experiments. Recent measurements of the 



SIDIS by the HERMES [38] and COMPASS [39] collaborations have shown clear 
evidence for the existence of the T-odd Sivers functions. These data also allow the first 
determination [40] of the magnitude and flavor structure of the Sivers functions and the 
nucleon transversity distributions. 

The TMD and transversity parton distributions can also be probed in Drell-Yan ex- 
periments. As pointed out [41] long time ago, the double transverse spin asymmetry 
in polarized Drell-Yan, Att, is proportional to the product of transversity distributions, 
hi{xq)hi{xq). The single transverse spin asymmetry, A^, is sensitive to the Sivers func- 
tion [42], f^{x) of the polarized proton (beam or target). Even unpolarized Drell-Yan 
experiments can be used to probe the TMD distribution function, since the cos2^ az- 
imuthal angular dependence is proportional to the product of two Boer-Mulders func- 
tions [43], hf {xi)hi {x2) . A unique feature of the Drell-Yan process is that, unlike the 
SIDIS, no fragmentation functions are involved. Therefore, the Drell-Yan process pro- 
vides an entirely independent technique for measuring the TMD functions. Furthermore, 
the proton-induced Drell-Yan process is sensitive to the sea-quark TMDs and can lead to 
flavor separation of TMDs when combined with the SIDIS data. Finally, the intriguing 
prediction [44] that the T-odd TMDs extracted from DIS will have a sign-change for the 
Drell-Yan process remains to be tested experimentally. 

No polarized Drell-Yan experiments have yet been performed. However, some infor- 
mation on the Boer-Mulders functions have been extracted recently from the azimuthal 
angular distributions in the unpolarized Drell-Yan process. The general expression for 
the Drell-Yan angular distribution is [45] 

da - 7 V 9 

— oc H-Acos 0-|-jUsin20cos0-|--sin^0cos2(^, (4) 
ail 2 

where 6 and are the polar and azimuthal decay angle of the in the dilepton rest 
frame. Boer showed that the cos 20 term is proportional to the convolution of the quark 
and antiquark Boer-Mulders functions in the projectile and target [46]. This can be 
understood by noting that the Drell-Yan cross section depends on the transverse spins 
of the annihilating quark and antiquark. Therefore, a correlation between the transverse 
spin and the transverse momentum of the quark, as represented by the Boer-Mulders 
function, would lead to a preferred transverse momentum direction. 

Pronounced cos2<^ dependences were indeed observed in the NAIO [47] and 
E615 [48] pion-induced Drell-Yan experiments, and attributed to the Boer-Mulders 
function. The first measurement of the cos 20 dependence of the proton-induced Drell- 
Yan process was recently reported for p + p and p + d interactions at 800 GeV/c [49]. 
In contrast to pion-induced Drell-Yan, significantly smaller (but non-zero) cos20 az- 
imuthal angular dependence was observed in the p + p and p + d reactions. While 
the pion-induced Drell-Yan process is dominated by annihilation between a valence 
antiquark in the pion and a valence quark in the nucleon, the proton-induced Drell-Yan 
process involves a valence quark in the proton annihilating with a sea antiquark in the 
nucleon. Therefore, the p + p and p + d results suggest [49, 50] that the Boer-Mulders 
functions for sea antiquarks are significantly smaller than those for valence quarks. 



FUTURE PROSPECTS 



Future fixed- target dimuon experiments have been proposed at the 120 GeV Fermilab 
Main Injector and the 50 GeV J-PARC faciUties. As discussed earUer, the Fermilab E906 
experiment will extend the d/ u asymmetry measurement to larger x region. Another goal 
of this experiment is to determine the antiquark distributions in nuclei at large x using 
nuclear targets. New information on the quark energy loss in nuclei is also expected. As 
discussed earlier, an advantage of lower beam energies is that a much more sensitive 
study of the partonic energy loss in nuclei could be carried out [21]. 

With the possibility to accelerate polarized proton beams at J-PARC [20], the spin 
structure of the proton can also be investigated with the proposed dimuon experiments. 
In particular, polarized Drell-Yan process with polarized beam and/or polarized target at 
J-PARC would allow a unique program on spin physics complementary to polarized 
DIS experiments and the RHIC-Spin programs. Specific physics topics include the 
measurements of T-odd Boer-Mulders distribution function in unpolarized Drell-Yan, 
the extraction of T-odd Sivers distribution functions in singly transversely polarized 
Drell-Yan, the helicity distribution of antiqaurks in doubly longitudinally polarized 
Drell-Yan, and the transversity distribution in doubly transversely polarized Drell-Yan. 
It is worth noting that polarized Drell-Yan is one of the major physics program at the 
GSI Polarized Antiproton Experiment (PAX). The COMPASS experiment at CERN will 
also measure -induced Drell-Yan on transversely polarized targets to extract Sivers 
functions from the Single-Spin-Asymmetry. Two recent Letter-of-Intent have also been 
submitted to RHIC for measuring polarized Drell-Yan at RHIC-Spin, one in collider 
mode and the other utiUzing internal targets. It is clear that the Drell-Yan process will 
continue to provide unique information on parton distributions in the forseeable future. 
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